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The electronic and magnetic properties of single 3d transition-metal(TM) atom (V, Cr, Mn, 
Fe, Co, and Ni) adsorbed graphdiyne (GDY) and graphyne (GY) are systematically studied 
using first-principles calculations within the density functional framework. We find that the 
adsorption of TM atom not only efficiently modulates the electronic structures of GDY/GY 
system, but also introduces excellent magnetic properties, such as half-metal and spin-select 
half-semiconductor. Such modulation originates from the charge transfer between TM adatom and 
the GDY/GY sheet as well as the electron redistribution of the TM intra-atomic s, p, and d or- 
bitals. Our results indicate that the TM adsorbed GDY/GY are excellent candidates for spintronics. 



PACS numbers: 73.20.Hb, 71.20.-b, 73.20.At 

I. INTRODUCTION 

Graphene ?i a single atomic layer of carbon sheet, 
has attracted tremendous attention because of its 
unique properties such as high room temperature 
carrier mobility-^ and quantum Hall effect Graphyne 
(GY), a two dimensional carbon allotrope with the 
same symmetry as graphene predicted to have a high 
possibility of synthesis by Baughman et al. in 1987 6 , 
is made up of hexagonal carbon rings and acetylene 
linkages. 7 So far, considerable experimental efforts 
have devoted to synthesize such unusual materials 
and numerous of related polymer structures have been 
prepared successfully £~— Graphdiyne (GDY), belonged 
to the family of GY, is also mesh structure consisting 
of carbon hexagonal rings and butadiyne linkages*^— 
Most intriguingly, (GDY) has been successfully grown on 
the surface of copper via a cross-coupling reaction using 
hexaethynylbenzene j 15 i 16 Consequently, the preparation 
of GY can be expected because theoretical work has 
shown that it is energetically more stable than GDY. 
In comparison with graphene, GY and GDY possess 
many interesting properties because they posses both 
sp-hybridized and sp 2 -hybridized carbon atoms and 
the natural 'holes'. For example, GY is a promising 
candidate for the anode material in lithium battery 
applications because of its higher lithium mobility and 
higher lithium storage capacity than graphite^ On 
the other hand, the Ca decorated GY is predicted an 
excellent candidate for hydrogen storage because of its 
high capacity hydrogen storage and effective prevention 
for the formation of Ca cluster J£ 

Analogous to the graphene, GDY/GY possess small 
carrier effective masses ; 12 : 19 high carrier mobilities^ 
In particular, based on density functional theory 
(DFT) calculations the GY and GDY are found to be 
semiconductors with direct band gap of 0.52eV and 
0.53eV, respectively^ 2 - which makes them advanced 
in preparation of effective room temperature field- 



effect transistors(FET). To expand the application of 
GDY/GY, modulation of their electronic properties is 
essential. Recent studies have indicated that similar to 
that of graphene, GDY/GY nanoribbons can effectively 
tailor the electronic properties of GDY/GY i 13 i 20 Inter- 
estingly, new families of carbon nanotubes based on 
graphyne motifs with even richer variation in electronic 
properties than do ordinary single wall nanotubes 
were predicate by Coluci et a l 21 ' 22 . Doping is another 
effective approach to modulate the electronic properties 
of GDY/GY. For example, GY with adsorption of 
Ca represents spin-polarized electronic property with 
a small localized magnetic moment(^0.25/iB).— It 
has been proved that the surface adsorption of TM 
atom profoundly influences the electronic properties 
of graphene and induces spin polarization 2 ^— in the 
system, which makes is a vital candidate materials for 
spintronics j 29 ' 30 However, to our knowledge, the effects 
of adsorption of TM atom on the surface of GY/GDY 
has not been reported as yet. It is significant to study 
such effect for their promising applications. To this end, 
we systematically investigate the electronic and magnetic 
properties of GDY and GY with single TM (V, Cr, Mn, 
Fe, Co and Ni) atom adsorption (TM-GDT/GY) using 
the first-principles calculations within the framework of 
the density functional theory (DFT). We find that the 
adsorption of TM atom not only efficiently modulates 
the electronic structures of GDY/GY system, but 
also introduces excellent magnetic properties, such as 
half-metal and spin-select half-semiconductor. 



II. METHOD AND CALCULATION DETAILS 

A super-cell of GY/GDY with TM adatom is adopted 
to be the calculated model as shown in Fig[IJ where a 
vacuum space of 13A perpendiculars to the GY/GDY 
plane is chosen to avoid the interactions between the 
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FIG. 1: Diagrams of TM-GDY/GY. (a) and (b) are the struc- 
tures of pure GDY and GY, respectively, (c) and (d) illustrate 
the most stable structures of TM-GDY and TM-GY, respec- 
tively. Possible adsorption sites are marked in (a) and (b), 
where A; ,Bi and Hi denote top, bridge and hollow sites, re- 
spectively. Blue and pink balls represent the carbon and TM 
atoms, respectively. 



neighboring images. The first-principles calculations are 
performed using the Vienna ab initio simulation pack- 
age (VASP)£Li22 within spin-polarized DFT, 3 ^ 3 ^ As re- 
ported by Wehling et ali^ 5 . that the electronic configura- 
tion, the adsorption geometries, and the magnetic state 
of the 3d-TM adatom on graphene are very sensitive to 
the effects of local Coulomb interactions U in the TM d 
orbitals. We employed a generxalized gradient approxi- 
mation functional (GGA)2& as well as GGA + U to de- 
scribe the exchange-correlation functional. We use U = 4 
and J=0.9 eV to study the dependence of the electronic 
structure on U, which is a typical value for 3d TM i 37 i 38 
The projected augmented wave (PAW) 3 ^^, potentials are 
adopted to describe the interaction between nucleus and 
valence electrons. A plane- wave basis set with the kinetic 
energy cutoff of 400 eV is employed. The Brillouin zone 
(BZ) is sampled using 5x5x1 and 9x9x1 Gamma- 
centered Monkhorst-Pack grids for the calculations of re- 
laxation and electronic structures, respectively. All sys- 
tems are fully relaxed without any symmetric constrains. 
The criteria of energy and atom force convergence are 
set to 10~ 5 eV/unitcell and 0.01 eV/A, respectively. The 
dipole correctional is considered to deal with the impact 
of variety of potential distribution introduced by TM ad- 
sorption. All calculation parameters in present work are 
systematically optimized. 

To evaluate the most favorable adsorption configura- 
tion, adsorption energy E a is adopted which is defined 



E a — E'_ 
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where ExM+sheet denotes the spin-polarized total energy 
of TM adsorbed GY/GDY sheet, E sheet is the total en- 
ergy of the isolated pure GY/GDY sheet, and Etm is 
the spin-polarized total energy of the corresponding free 
TM atom. 



To analyze spin polarization of the systems, we define 
the spin polarization P(Erf /homo)) a t the Fermi level 
or at the highest occupied molecular orbital (HOMO) as: 



P(E(f/HOMO)) = 



D{E(f/HOMO)t) — D{E(f/HOMO)l) 
E(E(f/HOMO)t) + D{E{f/HOMO)]) 



where D(E {f/HOMon ) and D(E {f /HOMO)i ) represent 
the DOS of the majority spin and minority spin at the 
Fermi level or the HOMO state, respectively. 

The charge transfer between the TM adatom and the 
GDY/GY sheet are evaluated using the atomic basin 
charge based on the atoms in molecular (AIM) theory^. 
The AIM theory gives a unique definition of partition- 
ing space into atomic basins by simply knowing the total 
charge density. The atomic basin is defined as a region 
of space traversed by trajectories of the electron den- 
sity gradient terminating at a given nucleus (attractor) 
and enclosed inside a zero charge-density flux surface 
^ p(r) ■ ~rt(r) = 0. Here ~rt{f) is the unit vector of the 
surface. Thus, space is divided into regions called atomic 
basin by the surfaces that run through the points where 
^ p(r) = 0. By integrating the electronic density within 
the basin, the charge of each atom can be estimated. 

To have a better understanding of the magnetism 
of the TM-GDY/GY system, we use the map of spin- 
polarized charge density(SCD) for all investigated sys- 
tems. The SCD are defined as: 

p{r) = p f (r) - p I (r) 

where p j" (r) and p 4- (r) are the spin up and the spin 
down charge density of TM-GDY/GY system, respec- 
tively. 



III. RESULTS AND DISCUSSIONS 

A. Adsorption geometry and energy 

Before investigating the TM adsorbed systems, we 
firstly optimize the pristine GDY and GY. The opti- 
mized lattice constants of GDY and GY are 9.50A and 
6.99A, respectively, which is in good agreement with pre- 
vious reports^—. In present work, six 3d TMs (V, Cr, 
Mn, Fe, Co and Ni) are selected to adsorb on the GDY 
and GY sheet. Considering the symmetry of GDY and 
GY, all possible initial adsorption sites are considered 
for each types of TM adatom, as shown in FigUJ (a) and 
(b). We find that the most favorite adsorption site is 
D x for GDY and D 2 for GY for all the six 3d TMs, as 
shown in FigfT] (c) and (d). All the TM adatoms dis- 
turb the lattice configuration of GDY/GY to some ex- 
tent. Figure [2] (a) and (b) give the TM adsorption ener- 
gies and structure parameters for TM-GDY/GY, respec- 
tively. From Figf2](a), one can see that the considerable 
negative adsorption energies for all the 3d TM adatoms 
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FIG. 2: (a) is the adsorption energy of TM-GDY/GY sys- 
tems, (b) is the distance between TM adatom and its nearest 
carbon atoms on the GDY/GY sheet and height between the 
adatom and GDY/GY sheet. Charge transfer (T(e)) from the 
TM atom to GDY/GY sheet based on AIM theory is shown in 
(c), and (d) is the magnetic moment (M) and spin-polarized 
energy AE 3pin . 



on GDY/GY indicating their chemisorption characteris- 
tics. The largest adsorption energy is -1.17eV (-1.39eV) 
when the Mn adsorbs on the GDY(GY) sheet; the mini- 
mum one is -2.87eV (-3.95eV) when the Ni adsorbs on 
the GDY (GY) sheet. From FigE] (b), we find that 
as for TM-GDY the maximum (minimum) distance be- 
tween the adatom and the nearest carbon atom (cLac) is 
2.282A (1.862A) for Cr (Co); and the maximum (mini- 
mum) height between adatom and GDY sheet is 0.135A 
(0.025A) for Fe (Ni). As for GY, the maximum (mini- 
mum) (Lac is 2.277A (1.948A) for Cr (Ni), which is sim- 
ilar to the case of GDY; whereas the maximum (mini- 
mum) height is 0.895A (0.043A) for V (Co). According 
to previous reports^—, the most favorable adsorption 
site for TM on graphene can be top, bridge and hollow 
site depending on the balance the contributions to chem- 
ical bonding from the s against those from the d orbitals 
of TMs2£, however, the most favorable adsorption site 
for TM on GDY/GY is all in the alkyne ring. Moreover, 
the adsorption energies for the TM-GDY/GY are much 
smaller than those of graphene, which indicates that the 
adsorption of TM atoms on the GDY/GY is more sta- 
ble than the corresponding TM adsorption on graphene. 
Such phenomena originate from the different hybridiza- 
tion of carbon atoms between graphene and GDY/GY. 
There are only sp 2 carbon atoms in the graphene but sp 
and sp 2 carbon atoms in the GDY/GY. The sp carbon 
atoms play a key role in the absorption process of TM 
atoms on the GDY/GY. Besides the sp 2 hybridization as 
in the case of graphene, which only exists out-of-plane 
p z ir/ir* (here z is the direction perpendicular to the car- 
bon plane), there are in-plane p x -p y ir/ir* in GDY or GY 
for sp hybridization. Such hybridization nature of GDY 
or GY enables the 7r/7r* to rotate toward any direction 
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FIG. 3: Spin-polarized charge density (SCD) distribution for 

V (a) , Cr (b), Mn (c), Fe (d) and Co (e) adsorbed GDY; 

V (f), Cr (g), Mn (h), Fe (i) and Co (j) adsorbed GY. The 
isosurface charge density is set to 0.004e/A 3 . 



perpendicular to the line of -C=C-, thus making it possi- 
ble for the 7r/7T* states from the -C=C- bonds at a given 
alkyne ring point toward the TM atom. As a result, 
the valance electrons of the TM atom couples with not 
only the p z but also the p x and p y of carbon atoms in 
the -C=C-, which produces the TM atoms adsorb in the 
alkyne ring with more favorable energy. Moreover, the 
heights between TM adatom and the GDY/GY sheet are 
all smaller than that of graphene sheet, which originates 
from both the bigger size of acetylenic and butadiynic 
ring and the rotation of 7r/7r* states of -C=C-. 



B. Magnetic and Electronic structures of 
TM-GDY/GY 

Based on the most stable structures of TM-GDY/GY, 
we systematically investigate their magnetic and elec- 
tronic properties. Distinct spin polarization is found for 
the GDY/GY with adsorption of V, Cr, Mn, Fe and Co. 
The corresponding magnetic moment is given in Figf2](d) 
and TablU The magnetism of TM-GDY/GY is mainly 
contributed by the TM adatoms. The sequence of mag- 
netic moment of TM-GDY/GY is Cr>Mn>V>Fe>Co 
for TM-GDY and Cr>Mn>Fe>V=Co for TM-GY, 
respectively. In comparison with the free-standing states 
of TM atoms, the magnetic moments of V, Cr, Mn, 
Fe and Co atom are all reduced when they adsorb on 
both GDY and GY. The reduction is derived from the 
electron injection from TM to GDY/GY and electron re- 
distribution between the electron orbitals of TM induced 
by the strong coupling between TM and GDY/GY, 
which reduce the number of spin-polarized electrons 
of TM atom and produce the TM adatom in low-spin 
configuration, the details will be discussed below. In 
Fig. [2] (d) and Tab. |H we also show the spin-polarized 
energy (AE sp i n ) which is the energy difference between 
the non-magnetic state and the magnetic state. All the 
AE spin s are larger than 0.7 eV, indicating the stability 
of the spin-polarized states of GDY/GY with adsorption 
of this five TM atoms, and some of them are hoping 
to be high temperature FM systems. The sequence of 
AE sp i n is similar to that of magnetic moment of the 
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TABLE I: The calculated results of TM-GDY/GY. "M","HM", "HSM" and "SM" in the band structures (BS) column denote 
metal, half-metal, half-semiconductor and semiconductor, respectively. T e is the charge transfer from the TM to GDY and GY 
derived from AIM. "4s*/3d*/4p*" denotes the the valence electron configuration corresponding to the Bader charge analysis 
of the TM-GDY/GY. TM^/TM^ is the electron numbers of spin up and spin down for TM atomic basin. Magnetic moment 
of TM-GDY/GY system /j,t((j,b) and magnetic moment of TM atomic basin /i a (/x_s) in TM-GDY/GY system are also given. 
AE S pi n is the energy difference defined as AE sp i n =ENM-EM , where Enm and Em are the total energies for non-magnetic and 
magnetic states in one unit cell, respectively. 
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V-GDY 


HM 


1.391 


0.269/3.187/0.153 


3.219/0.390 
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Cr-GDY 
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1.027 


0.402/4.290/0.281 
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Mn-GDY 
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1.047 


0.234/5.096/0.265 
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Fe-GDY 
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0.877 


0.417/6.243/0.463 
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2.16 
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Co- GDY 


M 


0.696 


0.494/7.324/0.517 


4.369/3.935 
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0.43 
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Ni-GDY 
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0.552 


0.388/8.677/0.383 


4.724/4.724 











V-GY 


HM 


1.567 


0.348/2.558/0.527 


2.286/1.148 


1.00 


1.14 
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Cr-GY 
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1.143 


0.394/4.102/0.361 


4.427/0.429 


4.58 


4.00 


2.975 


Mn-GY 


HSM 


1.260 


0.304/4.911/0.525 


4.203/1.537 


3.00 


3.05 


1.811 


Fe-GY 
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1.103 


0.451/5.819/0.626 


4.540/2.357 


2.62 


2.18 


1.949 


Co-GY 


HM 


0.731 


0.607/6.871/0.819 


4.742/3.555 


1.00 


0.19 


0.910 


Ni-GY 


SM 


0.572 


0.418/8.026/0.984 


4.714/4.714 












systems as a whole. For example, the Cr-GDY possesses 
the largest AE spin corresponding to the largest magnetic 
moment. Spin-polarized charge density (SCD) as shown 
in Fig. [3] represents that the coupling between TM atom 
(expect Ni atom) and GDY or GY is antiferromagnetism 
and the magnetic moment of TM-GDY/GY is mainly 
contributed by the TM adatom. Moreover, one can see 
that much stronger SCD locates around the adatom for 
Cr and Mn adsorbed GDY (GY), consistent with their 
larger magnetic moment. 

The spin-resolved band structures of TM-GDY and 
TM-GY are shown in Fig. 0] and Fig. [SJ respectively. 
Our results indicate that the pristine GDY and GY are 
direct band gap semiconductors with band gap of 0.49eV 
located at the T point and 0.48 eV located at the M 
point, respectively. The results are in good agreement 
with previous worki£. The adsorption of V, Cr, Mn, Fe 
and Co induces spin splitting of the band structures and 
produces the magnetism for TM-GDY/GY. However, 
Ni-GDY/GY shows no spin polarization and magnetism, 
which is consistent with the case of TM-adsorbed 
graphen o 23 ' 27 . Interestingly, versatile electronic charac- 
teristics are achieved in TM-GDY/GY. We summarize 
the magnetic and electronic properties of TM-GDY/GY 
as follows: 

(i) Half-metal: V-GDY, V-GY and Co-GY behave as 
half-metal which can be used as spin filter. Their band 
structures are given in Fig. [4] (a), Fig. [5J (a) and Fig. 
[5j (e), respectively. The V-GDY behaves as metal for 
majority spin channel and a semiconductor with an 
indirect band gap of 0.16eV for majority spin channel. 
The V-GY is metallic for majority spin channel and 
a semiconductor with an indirect band gap of 0.478 
eV for minority spin channel. However, the Co-GY is 
a semiconductor with a direct band gap of 0.426 eV 
for majority spin channel and behaves as a metal for 



minority spin channel. 

(ii) Spin-select half-semiconductor: The band structures 
of Mn-GY as shown in Fig. [5j (c) indicate its spin-select 
half-semiconductor with 100% spin polarization at 
HOMO state, where both the maximum of the valence 
band (VBM) and minimum of the condition band(CBM) 
belong to the same spin channel. The majority spin 
channel of Mn-GY is direct band gap semiconductor 
with band gap of 0.602eV at the M point. However, the 
minority spin channel of Mn-GY is indirect band gap 
semiconductor with band gap of 0.20eV at the K point, 
x (iii) Semiconductor-metal transition: The GDY and 
GY are spin-degenerated direct-gap semiconductors. 
However, when GDY is adsorbed by Cr, Mn, Fe and 
Co, as shown in Fig0] (b), (c), (d) and (e), and GY is 
adsorbed by Cr and Fe, as shown in FigJS] (b) and (c), 
the systems realize semiconductor-metal transition. All 
these systems are spin polarized metals. In particular, 
the spin polarization of Co-GY is up to 90%. 
(iv) Narrow-gap semiconductor: Although Ni does not 
induce spin polarization in GDY and GY, as shown in 
FigHl(f) and FigJS] (f), it narrows the band gaps of GDY 
and GY to 0.24eV and 0.40eV, respectively. The k point 
in reciprocal space of VBM and the CBM is the same as 
the pristine GDY/GY sheet. 

All the above results indicate that the adsorption 
of TM on GDY/GY is an effective modulation ap- 
proach on their magnetic and electronic properties. 
Such modulation mainly originates from the electron 
transfer between TM adatom and the GDY/GY, as 
well as between the TM intra-atomic orbitals induced 
by the strong coupling between TM and GDY/GY. To 
quantitatively investigate the modulation mechanism, 
we calculate the electron transfer for each adatom on 
the basis of Bader charge analysis.— The electrons 
occupied 4s, 3d, and Ap orbital of TM are obtained by 
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FIG. 4: Energy band structures of V-GDY (a), Cr-GDY (b), 
Mn-GDY (c), Fe-GDY (d), Co-GDY (e), and Ni-GDY (e). 
The solid (red) and dash (blue) lines represent the majority 
and minority spin channels, respectively. The Fermi level is 
set to zero. 
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FIG. 5: Energy band structures of V-GY (a), Cr-GY (b), 
Mn-GY (c), Fe-GY (d), Co-GY (e), and Ni-GY (e). The 
solid (red) and dash (blue) lines represent the majority and 
minority spin channels, respectively. The Fermi level is set to 
zero. 



separating total electrons in TM atomic basin based 
on the proportion of 4s, 3d, and 4p partial density of 
states (PDOS) (as shown in Fig|5] and FigJTJ) to total 
density of states (TDOS). Obvious charge transfer takes 
place from the TM adatom to GDY/GY sheet as shown 
in FigH (c) and TabUl For both GDY and GY, the 
maximum and minimum electron transfer occurs when 



FIG. 6: PDOS of V (a), Cr (b), Mn (c), Fe (d) , Co (e) and 
Ni (f) adsorbed GDY. The upper and lower panels are the 
PDOS of nearest neighbor carbon atom (CI) (as indicated in 
the[T](c)) of TM adatom and TM adatom, respectively. The 
Fermi level is set to zero. 



the systems are adsorbed by V and Ni, respectively. 
From TabU we find that the electrons not only transfer 
from TM to GDY/GY but also redistribute between the 
4s, 3d and 4v orbitals of TM adatom. The 4 s orbital 
of TM adatoms always lose electrons and the 4p orbital 
always get electrons. As for GDY system, the 3d orbital 
of V and Cr loses electrons, while the 3d orbital of other 
TMs gains electrons. However, as for GY system, only 
the 3d orbital of Co and Ni trivially gain electrons. 
Such charge transfer will affect the pairing of electrons 
of TM and further influence the magnetic properties 
of TM-GDY and TM-GY. The results in Tab ID show 
that the electrons occupied 4 s orbital of TM adatom 
for all TM-GDY/GY systems are all less than 0.6 e. 
Such orbital occupancy will produce the TM adatom on 
GDY/GY in the low-spin configurations as like the case 
of TM adsorbed graphene reported by Wehling et alj^. 

Such electron transfer is derived from the strong cou- 
pling between TM and GDY/GY which is determined 
by the bonding characteristics of GDY/GY. The p x and 
p y of sp hybridized carbon atoms in the GDY and GY 
contribute to both cr-bond and 7r-bond, making it easy 
to couple with the TM adatom. The projected density 
of states (PDOS) of TM-GDY/GY systems is shown in 
Fig|6] and Fig 13 to illuminate the coupling between the 
TM and the GDY/GY sheets. We only show the PDOS 
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FIG. 7: PDOS of V (a), Cr (b), Mn (c), Fe (d) , Co (e) and Ni 
(f) adsorbed GY. The upper and lower panels are the PDOS 
of nearest neighbor carbon atom (C2) (as indicated in the 
Fig[T](c)) of TM adatom and TM adatom, respectively. The 
Fermi level is set to zero. 



for one of the NN sp hybridized carbon atoms of the TM 
adatom in Fig [6] and Fig 13 because the TM adatoms only 
strongly bind to the nearest neighbor (NN) sp hybridized 
carbon atoms and the NN sp hybridized C atoms of the 
TM adatom in GDY or GY exhibit the similar results. 
Here we take Fe-GDY as example to elucidate the 
relationship between the coupling, electron transfer, and 
the magnetic properties of TM-GDY/GY. The PDOS of 
Fe adatom and CI for Fe-GDY are shown in Fig|6] (d). 
For the majority spin, the Fermi level passes through 
the conduction band due to the electron transfer from 
TM 3d to GDY p and the states around the Fermi 
level are mainly derived from CI p z states. Moreover, 
in the majority spin channel, the Fe d orbital with 
El (d yz ,d xz ) and E2 (d xy ,d x 2_ y 2) symmetry strongly 
couples with the p x and p z of CI in the energy ranges of 
-4.0^-3.4eV and -2.7~-1.0eV. For the minority spin, the 
states around the Fermi level are mainly derived form 
the coupling between Fe El and CI p z . One can see that 
there are considerable overlap for the Fe E2 and CI p x 
and p z in -3.5~-2.9eV and -1.6~-1.2eV energy windows. 
However, from the electron statistics based on the Bader 
Charge analysis^ as listed in Tab|H we find that though 
the Fe atom injects 0.877e electrons into GDY, the total 
electrons of Fe 3c? is still increased by 0.243e, which is 
derived from the electron promotion from Fe 4s to Fe 



3d. The intra-atomic electron transfer of Fe adatom is 
mainly determined by the coupling effect between the 
4p and 4s of Fe and CI p x and p z around -2.3eV, the 
coupling between the 4p, El, and E2 of Fe and CI p x 
and p z around -1.4eV in the majority spin channel, as 
well as the coupling between the Ap, El and E2 of Fe 
and p x and p z of CI around -3.0eV in the minority 
spin channel. Along with the coupling between Fe and 
GDY, the empty Fe Ap show 0.463e electron occupation. 
This is the result of electron exchange from GDY to 
the Fe adatom and the electron redistribution between 
intra-atomic orbitals from 4s and 3d to the Ap for Fe 
adatom. It can be found that the states of Fe Ap and 
CI p x couple with each other around -3.7eV, -3.0eV and 
-1.4eV in both majority and minority spin channels as 
shown in Fig[5] (d), which corresponds to the electron 
exchange from the CI p x to the Fe 4v states. Such 
electron exchange between Fe and GDY as well as the 
intra-atomic oribtal electron redistribution consequently 
affect the pairing of electrons of Fe and further determine 
the magnetic property of Fe-GDY. As a result, the Fe 
adatom behaves a low-spin configuration with 2.16/is 
magnetic moment in Fe-GDY system. The same physical 
mechanism results in the low-spin configurations of TMs 
in the TM-GDY/GY when GDY/GY is adsorbed by V, 
Cr, Mn, Fe, Co and Ni. 

Here, we summarize the electronic coupling character- 
istics near the Fermi level for TM-GDY/GY in details 
as follows: 

(i) Half-metal: The PDOS of half-metals V-GDY, V-GY 
and Co-GY are shown in Fig[5] (a), FigJTJ (a) and FigJT] 
(e), respectively. The V-GDY is metallic with the Fermi 
level passing through the conduction band determined 
by the states of p z of CI in the minority spin channel. 
As for the majority spin channel, the VBM is derived 
from the coupling of V El and CI p z , as shown in 
Fig|6] (d). FigJT] (a) indicates that the Fermi level of 
V-GY passes through the valence band in the majority 
spin channel where the VBM derive from the coupling 
between p x and p z of C2 and Al {d z -z) of V. However, in 
the minority spin, the VBM derives from the coupling 
between p x and p z of C2 and E2 of V and the CBM 
derives from the coupling between p x ,p z of C2 and El 
of V. As shown in Fig{7] (e), the Co-GY is a half-metal 
with the Fermi level passing through the valence band 
in minority spin channel which is mainly contributed by 
the coupling between the p z of C2 and El of Co. 

(ii) Spin-select half-semiconductor: Mn-GY is a spin- 
select half-semiconductor with 100% spin polarization 
at HOMO state. The PDOS of Mn-GY as shown in 
FigJT] (c) indicate that it is a semiconductor in both 
spin channels. In majority spin channel, both the VBM 
and CBM are derived from the coupling between p z and 
p y states of C2 and El states of Mn. In minority spin 
channel, the VBM mainly derives from the coupling 
between E2 states of Mn and p x states of C2, whereas 
the CBM derives from the coupling between C2 p z and 
Mn El. 
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(iii) Semiconductor-metal transition: The systems 
realize semiconductor-metal transition, when GDY is 
adsorbed by Cr, Mn, Fe and Co as shown in Fig|6] 
(b), (c), (d) and (e), and GY is adsorbed by Cr and 
Fe, as shown in FigJT] (b) and (d), respectively. All 
these systems are spin polarized metals with obvious 
spin polarization. The spin polarization of Cr-GDY, 
Mn-GDY, Fe-GDY, Co-GDY, Cr-GY and Fe-GY is 6%, 
18%, 20%, 90%, 53% and 20%, respectively. In all these 
systems, considerable p z of CI locate at the Fermi level 
in both spin channels. Obvious spin polarization of 
3d can be found in these systems. The p z of C locate 
around the Fermi level in the majority spin channel 
for Cr-GDY, Mn-GDY, Fe-GDY, Co-GDY and Fe-GY, 
respectively. The p z and El of Cr locate around the 
Fermi level in the majority spin channel for Cr-GY. 
The p z of C and El of Mn, Fe, Co locate around the 
Fermi level in the minority spin channel for Mn-GDY, 
Co-GDY, Fe-GDY and Fe-GY, respectively, whereas the 
p z of C locate around the Fermi level in the minority 
spin channel for Cr-GDY and Cr-GY. 

(iv) Narrow-gap semiconductor: Ni-GDY/GY narrows 
the band gap in comparison with that of pristine 
GDY/GY. As shown in FigJB] (f) and Figd (f), the 
majority and minority spin electronic states are degener- 
ated for both Ni-GDY and Ni-GY. And both the HOMO 
and LUMO derive from the p z states of sp hybridized 
NN carbon atoms and the El states of Ni. 



IV. CONCLUSION 

In summary, we have performed DFT calculations 
on the magnetic and electronic properties of single 



TM atom adsorbed GY/GDY sheet. We find that the 
adsorption of TM can effectively modulate the magnetic 
and electronic properties of GY/GDY sheet. V, Cr, 
Mn, Fe and Co adatoms induce the spin polarization 
and magnetism for GDY/GY. Most intriguingly, several 
interesting spintronics properties, such as half-metal 
and spin-select half-semiconductor, are achieved in the 
TM-GDY/GY. The modulation of the magnetic and 
electronic properties of GDY/GY by the adsorption of 
TM atoms originates from the coupling between TM 
adatom and GDY/GY and the electron redistribution 
between intra-atomic orbitals of TM. Our results indicate 
that TM-GDY/GY systems are potential candidates for 
spintronics. 
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